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ABSTRACT 

The supermassive black hole at the center of the Milky Way plays host to a massive, young cluster 
that may have formed in one of the most inhospitable environments in the Galaxy. We present new 
measurements of the global properties of this cluster, including the initial mass function (IMF), age, 
and cluster mass. These results are based on Keck laser-guide-star adaptive optics observat ions used 
to ide ntify the young stars and measure their Kp-band luminosity function as presented in IDo et ah! 
(|2012t ). A Bayesian inference methodology is developed to simultaneously fit the global properties 
of the cluster utilizing the observations and extensive simulations of synthetic star clusters. We find 
that the slope of the mass function for this cluster is a = 1.7 ± 0.2, which is steeper than previously 
reported, but still flatter than the traditional Salpeter slope of 2.35. The age of the cluster is between 
2.5-5.8 Myr with 95% confidence, which is a younger age than typically adopted but consistent within 
the uncertainties of past measurements. The exact age of the cluster is difficult to determine since our 
results show two distinct age solutions (3.9 Myr and 2.8 Myr) due to model degeneracies in the relative 
number of Wolf-Rayet and OB stars. The total cluster mass is between 14,000 - 37,000 M Q above 1 
Mq and it is necessary to include multiple star systems in order to fit the observed luminosity function 
and the number of observed Wolf-Rayet stars. The new IMF slope measurement is now consistent 
with X-ray observations indicating a factor of 10 fewer X-ray emitting pre-main-sequence stars than 
expected when compared with a Salpeter IMF. The young cluster at the Galactic center is one of the 
few definitive examples of an IMF that deviates significantly from the near-universal IMFs found in 
the solar neighborhood. 

Subject headings: galaxy: center - stars: luminosity function, mass function - stars: massive - stars: 
evolution - methods: statistical - infrared: stars 



1. INTRODUCTION 

Young nuclear star clusters have now been found 
surrounding supermas sive black holes in a number of 
nearby galaxies /e.g. lLauer et al.l 119981 : iBender et al.1 
12001 iSeth et al.l 12006^ The best studied young 
nuclear cluster is at the center of our own Milky 
Way Galaxy, located only 8 kpc away and sur- 
rounding a supermassive black hole (SMBH) of mass 
M~ 4 x 10 6 Ma, (lEckart fe GenzeH Il997l: IGhez et all 
19981 [20001. J2003l I2005L 120081: IGenzel et al.l 120001: 
"IT 



Schodel et aU2002l:lSchodel et all2003HEisenhauer et al l 



20051 : iGillessen et all I^Oga) . The young stars are lo- 
cated within 1 pc of the SMBH a nd are thought to hav e 
formed as recently as 4-8 Myr ago (jPaumard et al.ll2~006f ). 
The origin of the young stars, in such close proximity to 
the supermassive black hole, is puzzling, given that the 
strong tidal forces in this region will shear apart typical 
molecular clouds before they can collapse to form stars. 
Thus, the young nuclear cluster (YNC) is a key labora- 
tory for understanding whether and how stars form under 
such extreme conditions. 
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The YNC at the Galactic center has several observed 
properties that may help to determine its origin. To 
date, more than 100 young stars have been spectro- 
scopically identified as OB supcrgiants, Wolf-Rayet 
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population appears to fall into three dynamical cate- 
gories; (1) a well-defined clockwise rotating disk ranging 
from 0.03 pc to at least 0.5 pc with moderate eccentric- 
ities, (2) an off-thc-disk population in a more isotropic 
distribution over the same distances and also with 
moderate eccentricities, and (3) ^10% cluster within 
0.03 pc of the black hole with high eccentricities of 
e = 0.8 (IGenzel et al Il2000t iLevin fc Beloborodov|[200l 
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These kinematic differences have led to some speculation 
that the young stars may have been formed in multiple 
episodes, although the outer two group s (#1 and #2) 
have consistent stellar populations (|Paumard et all 
120061 ). The inner group (#3) was most likely dynam- 
ically injected, since the SMBH's tidal forces are too 
strong to permit star formation at these close distances. 
Initial theories suggested that this inner population 
was substantially older, as long times were needed to 
scatter inward a large number of binary systems that 
can interact with the SMBH and leave behind a star 
on a hig hly eccentric orbit with a small sem i-major 
axis (e.g. iBrown et all 120071 iPerets et all 120071 ) . How- 
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ever, more recent inclusion of additional dynamical 
effects suggests that the stars in the central region can 
be brought in more efficiently than initial ly thought 
(jLockmann et al.l 120091: iMadigan et all 120111 ). Thus, it 
remains uncertain whether the inner group of young 
stars was originally part of the outer groups and born 
in the same star formation event. In the outer groups, 
there is some observatio nal support for a po ssible warp 
in the clockwise disk (|Bartko et al.l 120091 ). a second 
face-on disk made of counter-clockwise rotating stars 
(jGenzel et alJ l200l iBartko et all Eopl and for sub- 
clust e rs of stars both on and off the disk (jMaillard et all 
l200l ISchodel et afl 120051 ILu et al.ll2005l) . although the 
statistical significance of th ese results is still debated 
(|Lu et all [20091: lYeldall2012f ). The total surface density 
profile of th e young stars in the plane of the sky is 
E oc i?" 1 (|Do et al.l 120121 ) and substantially steeper 
in the disk plane 5]^W disk K R • Early analysis 
of the bright stars (K<13) gave an age for the cluster 
of 6 ± 2 Myr based on the presence of Wolf-Rayct 
WR) stars and the p roportions of WR stars to O stars 
Paumard et al.ll2006D . In this same work, the observed 
number of WR and O stars and their luminosity function 
suggested that the the total cluster mass is 10 4 M Q and 
the initial mass function (IMF, dN/dm oc m ) is 
top-he avy with a slope significantly flatter than iSalpeteil 
(|1955l a = 2.35). However, these results were limited 
by the lack of sensitivity to less massive, main-sequence 
stars (K>13, M<20 M ). Deeper spectroscopic studies 
initially showed that the luminosity function for the 
young stars appears to have a sharp turn-over at 
K=13.5, suggesting that the present-day mass function 
is extremely top-he avy with a slope of a = 0.45 ± 0.3 
(|Bartko et al.ll2010T ). A top-heavy IMF is also supported 
by the Chandra X-ray observations of the region, given 
that lower mass young stars should still be coronally 
active and emitting cop ious X-rays that aren't detected 
(jNavakshin fe Sunvaevll2005ft . 

Several possible models for the origin of the young stars 
have been proposed and the two presently supported by 
observations include (1) formatio n in situ in a massive 
self-g ravitating molecular disk (|Levin fc Beloborodovl 
120031 ) or (2) disruption of an in- falling mas sive cluster 
formed much further away (jGerhard l200l . In situ 
formation models invoke the build-up of a gas disk 
surrounding the SMBH that reached a critical mass 
(~ 10 4 M Q ) approximately 6 Myr ago such that local 
self-gravity within the disk became sufficient to overcome 
the tidal shear, the disk c ollapsed along the vertical 
direction and formed sta r s (iLevin fc Beloborodovl 120031: 
Kolvkhalov fc Svunvaev| 119801 iShlosman fc Beeelman 
19891: [Morris fc Serabvnl 119961 iSan dersI 119981 i Goodman 
20031 iNavakshin fc Cuadral 12005ft . If such a disk was 
built up slowly, the gas and resulting young stars 
would largely be on circular orbits, which does not 
appear to be the case observationally (|Lu et alJ 120091 
IBartko et al.l 120091 ). Modified in situ formation models 
include the rapid infall of a single massive molecular 
cloud or the collision of two in-falling clouds to pro- 
duce non-cir cular orbits and the on-disk and off-disk 
populations ([Alexander et alJ 120071 ICuadra et al.1 120081 
Sandersl 1199a iVollmer fc Duschll 120011: INavak shin et alJ 
20071 iHobbs fc Navakshinll2009D . The surface density of 



stars formed in situ in a slowly-built gas disk may be 
steep (E oc i? -2 ); however, the surface density resulting 
from a cloud-cloud collision is less well constrained. For 
in-falling cluster scenarios, the young stars are formed 
a few parsecs away in a massive young star cluster 
that spirals in via dynamical fric tion and disrupts in 
the central parsec (|Gerhardll200l . In order for such 
a cluster to migrate into the central parsec within 
only a few million years, the cluster must be very 
massive and centrally concentrated and perhaps even 
host a n intermediate mass black hole (IMBH) at its 
center (IKim fc Morrisll2003l iPortegies Zwart et al.ll2003l 
McMillan fc Portegies Zwartl 120031 iGiirkan fc Rasiol 
2005 1 IHansen fc Milosavlievicl [20031 IKim et al.l I2004ft~ 7 
Such a scenario would produce the on-disk population 
and a flatter surface density profile (E oc i? -075 ), with 
the off-disk populatio n resulting from subsequent dy- 
nami cal perturbations (jHaas et al.ll201ll : IBaruteau et alJ 
120111) . In situ formation models are favored based on 
the density profiles and available time scales; however, 
no single model completely explains all of the observed 
properties. One important prediction of current in situ 
formation models is that the initial mass function should 
be very top-heavy due to the extreme temperatures, 
pressures, densities, a nd ambient radiation fi e lds present 
in th e central parsec (INavakshin et al.l 120061 iNayakshinl 
l2006t lAlexander et all 120071 120081 ICuadra et all 12008ft . 
However, a number of model parameters are still very 
uncertain (e.g. gas infall rate, temperature, pressure, 
and cooling time) and a wide variety of initial mass 
functions are still possible. 

The initial mass function, age, and mass of the young 
nuclear star cluster at the Galactic center is of particular 
interest both for understanding the origin of the young 
stars and exploring star formation under extreme con- 
ditions. Past observations of a turn-over in the near- 
infrared luminosity function and the low total X-ray 
emission apparently support a top-heavy IMF. How- 
ever , new observati ons presented in our companion pa- 
per, iDo et all (|2012L Paper I), show a luminosity function 
that does not turn over, but continues to rise; warranting 
a new analysis of the mass function. 

In this work, we compare the near-infrared photometry 
of the spectroscopically identified young stars in Paper 
I to models of star clusters in order to determine the 
age, mass function, and total mass of the young nuclear 
star cluster. Details on the observational sample and 
measurements are presented in [J2] Synthetic clusters are 
generated using stellar evolution and atmosphere models 
to produce individu al st ars a s th ey would appear at the 
Galactic center in ij3.ll In ^13.21 we present a Bayesian 
inference method for fitting the observed data with syn- 
thetic clusters to determine the cluster's properties and 
uncertainties. Results in $4] show that the present-day 
mass function is slightly flatter than a Salpeter IMF 
and rule out a mass function as top-heavy as previously 
claimed. We also find a younger age and our modeling re- 
quires the inclusion of binaries and multiple star systems 
to adequately fit the luminosity function. 

2. OBSERVED DATA & SAMPLE 

Observat ions and data used in this paper are described 
in IPaper J . Our data sample includes all stars with a 
non-zero probability of being a young, early-type star 
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(pyng > 0). For each star, the measurements consist 
of a Kp magnitude corrected for differential extinction 
(K'aa)i a -^-P uncertainty (<tk'), the probability of youth 
{pyng), and an indicator for whether the star is a Wolf- 
Rayct star. Only stars with K' AA < 15.5 are included. 
The modeling presented in this work utilizes these indi- 
vidual measurements rather than a binned Kp luminosity 
function (KLF). However for illustrative purposes, Fig- 
ure [1] shows our binned KLF compared with a model 
KLF for a cluster with an age of 6 Myr, extinction of 
Aks = 2.7, distance of 8 kpc, and an IMF slope of 
a = 0.45, consiste nt with previously published best-fit 
cluster parameters (|Bartko et al.ll20iol ). The discrepancy 
between our observed KLF and this model motivates the 
work presented below as we attempt to derive new best- 
fit cluster parameters and uncertainties. 

In our sample, we have chosen to include young stars 
at all radii, including the central 0'.'8 region immediately 
around the supermassive black hole. There are some sug- 
gestions that the stars in this central region were not 
formed in the most recent starburst, but rather in an 
earlier event. The brightest early- type star in this re- 
gion, SO-2, is a B0-B2.5 V star with an age of less than 
15 Myr bas ed on its measured te mperature, gravity, and 
luminosity dMartins et al.|[20 08i). The other early- type 
stars in this region h ave spectral features co nsistent with 
main-sequence stars ([Eisenhauer et al.l2005f ) and the old- 
est age for a main-sequence star at Kp=15.5 (Ak s = 2.7, 
d=8 kpc) is less than 20 Myr. Given these age constraints 
and the large uncertainties in the theoretical scattering 
efficiencies, it remains plausible that the young stars in 
the central region of our sample are the same age as 
the star s at larger radii. Furthermore, as described in 
iPaper J . including or discarding the stars within a ra- 
dius of 0'.'8 does not significantly change the shape of the 
KLF. 

3. MODELING THE OBSERVED CLUSTER 

While the early-type stellar population in the central 
parsec of the Galaxy c an be divided into several different 
kinematic sub-groups dPaumard et al.|[2006t iGhez et al 
20081 



Lu et all 2009; Bartko et al. 2009; Gillcssen et al 



20091)). it is not yet clear whether these sub-groups have 
different origins, or whether all the young stars formed 
from the same starburst and were differentiated via sub- 
sequent dynamical processes. Therefore, we choose to 
analyze the entire population of young stars as a sin- 
gle starburst cluster and to estimate the cluster's mass 
function, age, and total mass. 

The traditional approach to deriving a mass function 
would be to construct a binned KLF from the observed 
Kp magnitudes, assume a distance and extinction, and 
use a mass-luminosity relation from models of stellar evo- 
lution and atmospheres to convert from observed magni- 
tudes to initial stellar masses. This approach has several 
shortcomings, including arbitrarily choosing bin sizes, 
neglecting uncertainties on stellar brightness, not ac- 
counting for multiple systems, and correlations in the 
clust e r age, mass, and IMF dMaiz Apellaniz fe Ubedal 
120051: iMaiz Apellanizl 120081 : iMai'z-ApellanizI 120091 ) . More 
recent and rigorous statistical methods, particularly 
those using Bayesian inference, provide flexibility and 
produce robust estimates of the most p robable cluster pa - 
rameters and their uncertainties (e.g. lAllen et al.ll2005l ). 
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Fig. 1. — The observed Kp luminosity function (KLF) for early- 
type stars in the Galactic center nuclear cluster (Paper I) compared 
with previously proposed model KLFs. The observed Kp magni- 
tudes (red) are corrected for differential extinction to a common 
extinction value of A^ s = 2.7. A model KLF is shown for pub- 
lished best-fitting cluster parameters (blue) with an IMF slope of 
a = 0.45, age of 6 Myr, and a distance of 8 kpc l IBartko et al.l 
120101) . For comparison, a model KLF is also shown for a similar 
cluster with a Salpeter IMF slope of a = 2.35 (green). 

We have developed a Bayesian inference methodology for 
deriving the physical properties of an observed young star 
cluster using simulations of synthetic clusters. The pro- 
cess of simulating a synthetic young cluster is described 
in H3.ll and the full Bayesian methodology, implemen- 
tation, and testing is described in N3.2[ H3.3I and H'SAl 
respectively. Readers wishing to skip the detailed de- 
scription of the modeling may proceed to 531 where the 
resulting best-fit cluster properties and uncertainties are 
presented. 

3.1. Synthetic Clusters 

Model young star clusters are produced assuming an 
instantaneous starburst at time £, at distance d, with an 
observed cluster mass M c \^ s . The initial mass function 
of the cluster is described as a single power-law given by 
the probability density function 
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where a is allowed to vary and has the value a = 2.35 
in the case of a typical IMF fo r the high-mas s end of a 
solar- neighborhood population (|Salpeterlll955[ ). In addi- 
tion to these free parameters, a number of fixed param- 
eters are also needed to fully model the cluster. First, 
the extinction is fixed to A^ s = 2.7 as we will com- 
pare to observations that have been corrected for differ- 
ential extinction to this value. The metallicity is fixed 
to Z=0.02, which is roughly solar. Observations sug- 
gest the young GC population is consistent with so- 
lar iron abundance; but may be as high as 2x solar 
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and seems enhanced in q-elements (iCarr eTaTl [20001: 
Ramirez et all 120001: iNaiarro et al.l 12004 iCunha et al l 
20071: iMartins et al] 120081 ) . The choice to keep metal- 



licity fixed to solar is also motivated by the lack of 
complete stellar evolution and atmosphere models at 
other metallicities and abundance ratios. Section [5] dis- 
cusses possible uncertainties associated with restricting 
our analysis to solar metallicity. The range of masses 
we consider is set as m m i n = 1 M and m max = 150 
Mq and is kept fixed throughout the analysis. This is 
justified as our observations do not include stars be- 
low 10 M Q due to sensitivity limits and stars above 150 
M Q have already exploded for clusters older than 3 Myr. 
Thus the currently observed stars have no power to con- 
strain the low or high mass cutoffs and a more expan- 
sive mass range would not impact our results. For each 
star drawn from the I MF, we must also conside r whether 
it has companions (ISagar fc Richtlerl flOQlL IKroupal 



ii uas co mpanions uoa.t'ai ca, ruc mieii uatm iivioiiua 
1995T: IGoodwin fe Kroupal 120051: iThies fe Kroupal I2007T 



Mafz- Apellanizl 120091: IWeidner et al.l I2009D . High mass 
stars in the Galactic disk are known to have a high multi- 
plicit y fraction, perhaps even 100% (jKobulnickv fc Frverl 
120071 ) and the multiplicity fraction and the companion 
star fraction, which is the mean number of companions 
per primary, are dependent on t he mass of the primar y 
star in the multiple system (e.g. lLafreniere et al.ll2008D . 
The Galactic center cluster is modeled allowing multi- 
ple systems with a mass-dependent multiplicity fraction 
(MF) and companion star fraction (CSF) . The functional 
form of how the MF and CSF vary with primary mass 
is still very uncertain even for solar-neighborhood young 
clusters. We empirically derive this relation by compil- 
ing MF and CSF measurements in young clusters from 
the literature and assuming that the MF and CSF fol- 
low a power law distribution (see Appendix [A] for the 
derivation). Each star drawn from the IMF with mass m 
is determined to be in a multiple system and assigned a 
number of companions based on 



MF(m) = 0.44 



CSF{m) = 0.48 
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always < 3. (4) 



The mass of the companion stars is drawn from a power 
law probability distribution on the mass ratio, q: 



P(l) 



1 + /? 



1 



(5) 



where /3 = —0.4 and the lowest allowed mass ratio 
is qi = 0.01. Our choice of j3 is consistent with 
observations of local massive stars in young OB as- 
sociations, although values ranging from f3 = 0.4 — 
—0.6 would also be within 68% confidence in t ervals 
allowed by observation^ JKobulnickv fc Frverl 120071 : 
iKiminki fc Kobulnickvll2012l ). Luminosity functions are 
hardly impacted by this choice of qi since lower mass 
ratios do not significantly increas e the total system lumi- 
nosity l|Kouwenhoven et al]|2009l ). Allowing for multiple 

6 The resulting companion star masses are substantially higher 
than randomly drawing companions from the IMF. 



systems will have the largest effect at the bright end of 
the luminosity function where magnitude bins that were 
un-populatcd by individual stars are filled in with mul- 
tiple star systems composed of massive, bright primaries 
and companions (|Maiz-Apellanizl 120091: IWeidner et"aH 
120091 ). In M3.4L simulated clusters are used to show that 
incorporating multiple systems is essential and that fail- 
ing to include them produces poor fits and incorrect ages 
and IMF slopes. For a cluster with a desired cluster mass 
of Md.obs, stars are simulated until the sum of all stellar 
masses, including companions, exceeds M c ; j0 f, s and then 
the last simulated star is thrown out if that would bring 
the synthes ized cluster mass closer to the desired M c ^ i, s 
(|Haasll20"l0h . A complete list of free and fixed parameters 
we use to produce synthetic star clusters is presented in 
Table [Q 

Conversion from stellar masses to synthetic photom- 
etry requires the use of stellar evolution and atmo- 
sphere models. At the suspected age of the young 
cluster (4—8 Myr), stars with masses in excess of 30- 
60 Mq are evolving off the main-sequence and stellar 
evolution models with rotation are found to best de- 
scribe populations of blue sup ergiants and Wolf-Rayct 
stars (|Mevnet fc Maederl [2003h . We use the grid of 
Geneva models with an initial rotation s peed of 300 
km/s provided bv lMevnet fc Maederl (|2003f ). which only 
extends down to 9 M Q. In these models, Wolf- 
R ayet stars are ident i fied u sing the criteria suggested 
in iMevnet fc Maederl ((20031 ): T efS > 10 4 K and the 
mass fraction of hydrogen at the surface is X s < 0.4. 
To extend to lower masses along the main sequence 
and also capture pre-main sequence evolution at younger 
ages, the Gene va grid is merged with the grid of mod- 
els provided by iSiess et all (|2000f) that extends up to 7 
M Q . For both grids, a suite of isochrones is generated 
at logarithmic ages of \ogt = 6.0 — 7.4 with steps of 
Alogi = 0.01 and each isochrone samples mass more 
finely than Am = 0.005 Mq . The two isochrones are 
then merged and interpolated over the gap between the 
two models (7-9 M Q ). An example of the merging pro- 
cess for logi = 6.8 is shown in Figure [5] The resulting 
isochrones provide a mapping from the simulated stellar 
masses to their luminosities (L), effective temperatures 
(T e ff) and surface gravities (g). 

The physical parameters for each star are converted 
into observable brightness at Kp-band using model stel- 
lar atmospheres. Again, there is no single set of at- 
mosphere models that spans the full range of effective 
temperatures and surface gravities. For hot stars with 
T e ff> 7000 K, atmosphere models bv lCastelli fc Kuruczl 
(|2003f ) are used for all non-Wolf-Rayet stars. For temper- 
atures of 4000 K <T e ff < 7 000 K , we use the NextGen 
models of lHauschildt et all (fl999h and for T eff < 4000 
K we use the up dated NextGen mo dels with improved 
AMES opacities (|Allard et al.ll2000D . The grids do not 
intersect exactly at the transition temperatures; however, 
synthetic broadband photometry differs by <2% in the 
K-band, which is less than our photometric precision, so 
no interpolation between the different atmosphere mod- 



7 IEkst rom et al. (2012]) contains an updated grid of Geneva mod- 
els with rotation extending down to lower masses and with rotation 
speeds other than 300 km/s; however, they were not available at 
the time of our analysis. 
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Fig. 2. — 6 Myr luminosity-effective temperature isochrones from 
various stellar evolution models used to construct and test a com- 
bined grid of stellar evolution models that spans ages from 1-25 
Myr and the full range of stellar masses. Geneva models without 
rotation spanning .8—120 Mq are shown in green as a baseline 
IjSchaller et al.iri992T ). The observed post- main-sequence evolution 
of massive stars at this age is best described with models that 
include rotation. The isochrone from Geneva models with rota- 
tion, shown in blue, differs significantly in luminosity. Models with 
rotation only exten d down to 9 Mq (blue star and green star). 
Models provided by Sicss ct al. (2000) arc used for main sequence 
and pre-main sequence models below 7 Mq and are shown in red 
(7 Mq point also shown as a red square and green square). The 
blue and red isochrones are merged and interpolated over the gap 
to construct final isochrones, shown in orange, used in our cluster 
modeling. 

els is done. All downloaded atmosphere models were 
re-formatted to work with the pysynphot python pack- 
age in order to perform interpolations within each grid 
and generate synthetic spectra at specific temperatures, 
mctallicitics, and gravities. The synthetic spectra are 
flux calibrated to what would be observed at Earth with- 
out extinction by multiplying by (R/d) 2 , where R is the 
radius of the star and d is the distance. Then the syn- 
thetic spectra are reddened using the Gal actic Center ex- 
tinction law from iNishivama et al.l (|2009() and convolved 
with transmission profiles for a typical Mauna Kea at- 
mosphere and for the Kp filter in NIRC2. The final syn- 
thetic spectra are integrated using pysynphot to produce 
broad-band Kp photometry as would be seen from the 
telescope. An example mass-Kp magnitude relation is 
shown in Figure [5] (left panel) and would predict that 
SO-2, with a Kp= 14.39 after correcting to Ak s = 2.7, 
would have a mass of ^18 M Q , which is consistent with 
previous estimates based on the effective temperature 
and surface gravity der ived from detail spectral analy- 
sis (jMartins et al.H2008l) . 

Wolf-Rayct (WR) stars require special consideration 
as their atmospheres are windy and clumpy. Mod- 
els of stellar evolution for WR stars can predict when 
stars enter and exit the WR phase, surface abun- 
dances, and the luminosities and temperatures at the 
base of the atmosphere. However, this is insufficient to 
uniquely define what will be observed due to differences 
in wind velocities, mass loss rates, and clumping factors 



(e.g. lHamann et al.l 12008ft . Put another way, there is 
no unique mass-magnitude relationship for Wolf-Rayct 
stars. Model atmospheres do exist for Wolf-Rayet stars; 
but they require observations of the wind velocities and 
mass-loss rates and still do not uniquely predict the inte- 
rior properties of the star, such as t emperature and gray - 
ity at the base of the atmosphere ([Martins et al.l 120071) . 
Therefore, the observed magnitudes of the Wolf-Rayet 
stars cannot be used in a straightforward manner and we 
exclude them from the analysis of the luminosity func- 
tion. The number of WR stars, N\yn relative to the 
number of OB stars is still a powerful constraint on the 
age and mass of the cluster and we incorporate Nvi/.r into 
our Bayesian analysis. 

3.2. Bayesian Methodology 

A Bayesian inference approach is used to determine the 
physical properties of the YNC at the Galactic center. 
This methodology is generally applicable to the analy- 
sis of other star clusters. The advantages of Bayesian 
in ference for a stronomical applications are highlighted 
in IPressj (119971) and have b een widely adopted in cos- 
mology (jHobson et al.l [2010h . Bayesian techniques are 
beginning to be used more widely in the study of all 
types of star clusters (ISelman et al.l 119991: lAllen et all 
20^IConverse fc Stahlerll2008t iDe Gennaro et alj|2009t 



van Dvk et alll2009HEspinoza et aLll2009HRizzuto et all 
20111) . The analysis presented here differs from previ- 



ous work in several ways. First, single-band photometric 
information is combined with some spectroscopic infor- 
mation that identify which stars are Wolf-Rayct stars. In 
our case, only single band photometry is necessary since 
extinctio n maps are available from past multi-band pho- 
tometry {Schodel et al. 2010). Second, the initial mass 
function is assumed to have the functional form of a trun- 
cated power-law, although other functional forms can 
easily be inserted. This differs from some Bayesian treat- 
ments that attempt to derive the masses of the individual 
stars allowing for infinite variation in the shape of the 
mass function (IConverse fc Stahlerll200l I van Dvk et all 
[2001 IDe Gennaro et alll2009D . Third, we allow for mass- 
dependent multiplicity fractions and high-order multi- 
ple systems, which are important at the high mass end 
of the mass function. Fourth, we include all candidate 
young stars and weight them by their probability of 
youth in order to use all available information and ac- 
curately account for incompletenes^]. Finally, instead of 
using Markov Chain Monte Carlo techniques to perform 
the Bayesian inference, we find that be tter solutions re- 
sult from nested sampl ing techniques ([Feroz &: Hobsonl 
120081: iFeroz et all [2009?). This Bayesian approach pro- 
vides detailed probability distributions for the mass func- 
tion slope and other cluster parameters given limited ob- 
servations. It also avoids traditional biases from bin- 
ning and neglecting photometric errors and accounts for 
stochastic sampling of stellar masses. Our full method- 
ology is described below in more detail. 
We start at the beginning, with Bayes equation 



P (0\D) 



P (D\e) P (Q) 
p(D) 



(G) 



A significant number of the candidate young stars were 
spectral-ty ped manually a nd were assigned 100% probability of 
youth. Sec Do ct al. (2012) for a complete discussion. 
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TABLE 1 

Model Cluster Parameters 



Parameter Vary? Description Prior 





free 


llVli. 3lUpL UtLwccll ''"THin dll*-l ' ' l //7 (J X 


TTnifnrm from TD 1 fl 3 

\J llllOl 111 11 Olll o.ou 


log £ 


tree 


logiQ^ age of the young cluster in years ) 


Gaussian with logt — 6.78, 








from [6.20, 6.70] 


^cl,obs 


free 


initial mass of the young cluster 


Uniform from [10 3 , 10 5 ] M 


d 


free 


distance to the young cluster 


Gaussian with d = 8.096 kpc, 






kpc from [6.793, 9.510] 




fixed 


IMF minimum stellar mass 


1 M 




fixed 


IMF maximum stellar mass 


150 M Q 


z 


fixed 


mctallicity of the young cluster 


0.20 (roughly solar) 




fixed 


extinction to the young cluster 


2.7 


MF_amp 


fixed 


MF function's amplitude 


0.44 


MFJndcx 


fixed 


MF function's power-law index 


0.51 


CSF_amp 


fixed 


CSF function's amplitude 


0.50 


CSFJndex 


fixed 


CSF function's power-law index 


0.45 


CSF_max 


fixed 


maximum value for CSF 


3 




fixed 


power-law index for q distribution 


-0.4 



where D is the observed data, is the model, p(D\Q) is 
the likelihood of observing the data given some model, 
p(Q) captures the prior knowledge on the model param- 
eters, and p{D) is a normalizing factor known as the 
evidence. The likelihood function and the prior distribu- 
tions then give posterior probability distributions for the 
model parameters, p(Q\D). The model, 0, is defined by 
the set of free and fixed parameters given in Table [TJ 

The data, D, is constructed from the sample of can- 
didate young stars (p yn g > 0) brighter than Kp=15.5 
as described in $2] The WR stars are separated out and 
only the observed number of WR stars, N WR , is included 
in our analysis, since the initial masses of WR stars can- 
not be determined based on their luminosity alone. For 
all the other non-WR stars, measurements of the Kp 
brightness and error are used along with the probabil- 
ity that each star is young, p yny - Although wc observe 
N i, s non-WR stars, not all of them are young. We esti- 
mate the number of young OB stars, Nob, as the sum 
of the observed non-WR stars, weighted by their prob- 
ability of youth: Nob = Y^i=i° Pyng,i- The resulting 
measurements that make up the data are then 

D =[N WR ,k obs ,N B] (7) 

where k Q b s is the set of {ki, 07^ , p yn g,i} measurements 
for the N i, s non-Wolf Rayet stars observed. 

The likelihood function is composed of three indepen- 
dent probabilities: (1) the likelihood of observing the set 
of Kp magnitudes, {ki}, and uncertainties, {0^}, which 
effectively captures the shape of the KLF; (2) the like- 
lihood of observing Nob OB stars, which captures the 
normalization of the KLF; and (3) the likelihood of ob- 
serving Nwr Wolf-Rayct stars: 

p(D\Q) =p(k obs |8) -piNoBlO). -p(N WR \Q). (8) 

To derive the first term, we use an unbinncd approach 
such that the likelihood is the multiplication of the in- 
dividual stars' probabilities, weighted by the probability 
of youth, 



N„ 



P (k obs |0) = J] P (h\e)Py- 



(9) 



given by the probability distribution derived from syn- 
thetically "observing" a simulated young cluster. First, 
the cluster is simulated with model parameters describ- 
ing the age (i), IMF slope (a), cluster mass {M c ^ bs), 
and the cluster distance (d), 

p(k\Q) intrinsic = Simulated Yng Cluster^, a, M c i t0bs , d). 

(10) 

These free parameters along with their priors are given 
in Table [TJ The resulting Kp photometry from the simu- 
lated stars in the model cluster are binned finely (0.1 mag 
bins) to produce a probability distribution for the intrin- 
sic Kp luminosities. For small numbers of stars, stochas- 
tic sampling effects can lead to an inaccurate estimate of 
the probability distribution. Therefore, in order to max- 
imize accuracy and reduce the total number of clusters 
that need to be simulated, all model young clusters are 
simulated with a total mass of 10 6 M Q . The shape of 
the Kp luminosity function does not change with clus- 
ter mass and the number of Wolf-Rayet and non-Wolf- 
Rayet stars simply scales linearly with total cluster mass. 
The intrinsic probability distribution for Kp luminosities 
is multiplied by the imaging completeness curve, C(k), 
truncated at Kp=15.5 to match the data, and normalized 
to give the final probability distribution, 



observed 



intrinsic 

C(k) 

Jo P(k\0) intrinsic 

C(k) dk 



(ii) 



Finally, the measured photometric error is incorporated 
by modeling the observed magnitude as ki = k' + e where 
e is drawn from a normal distribution with zero mean and 
a standard deviation of . The resulting probability of 
observing ki for a given star is then given by 



J){k I observed 



1 



e "ii dk', 



2 ™l 



The probability for a star to have some Kp, p(ki\Q), is 



(12) 

which is the convolution of the model KLF and the er- 
ror probability distribution. The resulting probabilities 
for each star feed into Equation [5] to calculate the first 
term of the total likelihood in Equation [5] In the second 
and third terms of the likelihood, the number of OB and 
WR stars scales linearly with the mass of the young clus- 
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ter, which was verified with simulations. Thus the same 
simulated young clusters used to produce p(k\&) intrinsic 
also give the expected number of WR and OB stars after 
scaling linearly with cluster mass, 



(13) 



VlO 6 M 

Xo B (M cl , obs ) = K O b(1O 6 M ) (^gg) (14) 

where t^wR and ^ob are the expectation values (means) 
for the number of WR and non-WR stars, respectively. 
The likelihood of observing Nwr and Nob are then 
taken as a Poisson distributions 

P[N W r\Q) = ^— j (15) 

p(iv OB |e) = — j . (16) 

3.3. Sampling Posterior Probability Distributions with 
Multinest 

The resulting posterior probability distribution for 
the model parameters cannot be calculated analytically. 
Traditionally, Monte Carlo techniques are used to pro- 
duce a representative sample of points from the poste- 
rior probability distribution. The most commonly used 
method is the Markov Chain Monte Carlo (MCMC). 
However, MCMC methods may have difficulty converg- 
ing or fully mapping probability space when the prob- 
ability distributions that are being sampled are multi- 
modal or highly degenerate and c urved. An a lterna- 
tive method is nested sampling (jSkillingl I2004T ) such 
as in the pub licly available multi-modal nested sam- 
pling code, \MultiNest\$FeTOz fc Hobsonll2008tlFeroz et all 
120091) . This algorithm has been successfully demon- 
strated in cosmolog y, galaxy evolution, and gravitational 
wave problems (e.g. [Bridges et al.ll2009t iKilbinger et al.l 
[20101 [Martinez et al.1 1201 H IVeitch fc Vecchidl2010f ) and 



is less computationally expensive and more accurate, 
in some cases, than MCMC methods. We performed 
tests on simulated clusters using traditional MCMC tech- 
niques using bo th Metro polis-Hastings and Hit-and-Run 
step methods (PyMC package) and compared the over- 
all computation time and accuracy with the MultiNest 
code. MultiNest took 5-10 times less computation time 
and produced more accurate results in all our simulated 
clusters (see H3.4l for more details). Therefore, we utilized 
MultiNest in our analysis of the Young Nuclear Star clus- 
ter's IMF. 

MultiNest uses a fixed number of points per iteration 
to sample parameter space and calculate the Bayesian 
evidence at each point. In successive iterations, the 
same number of points are concentrated into smaller and 
smaller volumes centered around the most probable re- 
gions. This process continues until the evidence no longer 
changes by more than a specified tolerance value. We use 
1000 points and an evidence tolerance of 0.3 in order to 
sample parameter space well and run efficiently. We also 
enabled multi-modal searches as several simulated clus- 
ters could clearly be fit by several distinct sets of model 
parameters. 



3.4. Testing on Simulated Clusters 

To test the speed and accuracy of our bayesian method- 
ology and codes, clusters were simulated and syntheti- 
cally "observed" that had different ages, IMFs, and mul- 
tiplicity properties. The synthetic "observations" were 
then fit with the bayesian inference techniques described 
above and the resulting posterior probability distribu- 
tion functions were examined to see if they produced re- 
sults consistent with the input cluster parameters. All 
clusters were simulated at a distance of 8 kpc, an ex- 
tinction of Aks = 2.7, and with a cluster mass aimed 
at producing a similar number of OB stars to our ob- 
served data set. Photometric errors for simulated cluster 
stars are randomly drawn from the distribution of em- 
pirical errors from the GC observations. Tests were per- 
formed to (1) understand degeneracies in fitting a clus- 
ter's age, (2) determine the impact of multiplicity, and 
(3) explore whether we recover both top-heavy and nor- 
mal IMF slopes. 

3.4.1. Age 

A cluster's age is often one of the most difficult param- 
eters to constrain, especially at young ages and when 
only high-mass main-sequence and post-main-sequence 
stars are observed. The number of WR stars and the ra- 
tios of WR to OB stars can be a precise indicator of 
age; however, there are also cas es where degeneracie s 
produce several possible solutions (jPaumard 
To illustrate these degeneracies we simulate clusters at a 
range of ages from 1.5 - 10 Myr and examine the ratio 
of Nwr to No b ■ Note that No b includes all OB stars 
down to Kp=15.5, which is equivalent to Bl— B2 V. Fig- 
ure[3]shows that ^wr/^OB initially rises around 2 Myr, 
peaks at 3 Myr, and falls to a minimum around 6 Myr, 
which is set almost entirely by the rise and fall of Nwr- 
After 6 Myr, the ratio rises slightly again as Nwr re- 
mains relatively constant, but Nob begins to drop. If we 
observe a Nwr/Nob ratio of 0.125, then there are two 
possible ages: 2.9 Myr and 3.8 Myr, assuming an IMF 
slope of a — 1.7. This degeneracy is not entirely indepen- 
dent of the IMF slope (Figure El right panel). In order 
to explore how our bayesian methodology handles this 
age degeneracy, we produce ten different realizations of 
a cluster, all with the same input parameters (t=4 Myr, 
a = 1.7, Md.obs — 10 4 M , d= 8 kpc, single stars). Each 
cluster is then fit using the Bayesian inference methods 
described in ^3.21 The resulting probability distributions 
for the age are shown in Figure [4] and multiple age solu- 
tions are almost always found. We note that in these ID 
marginalized probability distributions, the input cluster 
age is not always the highest peak. However, the in- 
put age is always recovered within the 99% confidence 
interval and with much higher confidence if each peak 
is treated as an independent solution in a multi-modal 
distribution. Given the results of our testing, all further 
age results will be reported using multi-modal solutions, 
if necessary. 

3.4.2. Multiplicity 

The need to incorporate multiplicity when determin- 
ing a cluster's initial mass function is often discussed 
(IKrouDalll995t IGoodwin fc Kroupal [20051 IWeidner et al.1 
l2009f ) but rarely implemented due to incomplete knowl- 




Age [Myr] 



Fig. 3. — Ratio of the number of WR stars to the number of 
OB stars (Kp<15.5) as a function of cluster age. This quantity, 
along with the magnitude of the brightest observed OB stars, pro- 
vides the strongest constraint on the cluster age; however, it shows 
degeneracies that often result in multiple solutions for a given clus- 
ter. Left: An example cluster with an IMF slope of a = 1.7 will 
produce a specific Nwr/Nob ratio two to three times in the first 
10 Myr (blue solid). Right: Changing the IMF slope will shift the 
ages at which a specific J^wr/^OB ratio occurs. The observed 
Nwr/Nob ratio for the YNC intersects the a = 1.7 cluster at 
ages of ~3 and ~4 Myr. Uncertainties due to stochastic sampling 
(shaded) also suggest that 8-10 Myr may be a plausible solution; 
but, the presence of bright O giants and supergiants makes these 
larger ages unlikely. 

edge of multiplicity fractions, number of companions, 
companion mass functions, and how these all scale with 
primary mass (see Appendix [X] and references therein) . 
The presence of multiple systems influences the bright 
end of the KLF as illustrated in Figure [5] by populating 
the brightest magnitude bins (Kp<ll) with sources that 
would otherwise be much fainter single stars. We tested 
two possible scenarios to explore the impact of multiplic- 
ity in our analysis. First, a cluster was simulated with 
only single stars and analysis was done with and without 
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Fig. 4. — Test-clusters' age posterior probability distributions 
for 10 different realizations of a cluster with a = 1.7, t=4 Myr, 
d=8 kpc, M c i j, s = 10 4 Mq , and only single stars. The input 
age of 4 Myr is marked (vertical red) and the confidence interval 
that just includes 4 Myr is reported in the top right corner of each 
panel. The confidence interval is calculated by sorting the prob- 
ability distribution and integrating starting with the highest bins 
and proceeding until the bin including 4 Myr is reached. The input 
age always falls within the 99% confidence interval. If each peak is 
treated as an independent solution, the input age is recovered with 
even greater confidence. 

allowing multiple systems in the fit. Second, a cluster was 
simulated with multiple systems and, again, analysis was 
done with and without allowing multiple systems in the 
fit. Both clusters had an IMF slope of a = 2.35, an age of 
4 Myr, and a mass of 10 4 M , although similar tests were 
performed with different IMF slope and age combinations 
and the results arc robust. Figure [6] shows the results for 
the simulated cluster containing only single stars. Incor- 
rectly fitting this cluster with multiple systems produces 
extremely biased estimates for the cluster age and only 
slight biases in the IMF slope and cluster mass. Without 
prior information, it is difficult to discern that an inac- 
curate multiplicity model is being used. Figure [7] shows 
the results for the simulated cluster containing multiple 
systems. Incorrectly fitting this cluster with single stars 
produces extremely biased estimates for nearly all the 
parameters. Even without prior information, it is easy 
to discern that the single-star model is inaccurate since 
the distance tends toward unphysical low values. Multi- 
ple systems in clusters produce Kp magnitudes that are 
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Fig. 5. — Comparison between the Kp distribution for model 
clusters with only single stars (green) and with multiple systems 
(red). Both clusters have the same basic parameters: t =4 Myr, 
d =8 kpc, A Ks =2.7, and M c i y = 10 6 M . Left: The mass- 
magnitude relation for two simulated clusters showing that multi- 
ple systems with high mass primaries are typically 1 mag brighter 
than single stars. Right: The binned Kp luminosity function for 
both clusters. The cluster with multiple systems extends to higher 
Kp magnitudes than does the cluster with only single stars. We 
also note the bump at Kp~18 is the pre-main-sequence turn-on. 

brighter than is allowed in the model clusters with only 
single stars at 8 kpc. Thus the only way to produce such 
bright systems is to move the cluster closer in distance. 
Shifting the cluster to younger ages is also a possibility; 
but such young clusters do not give the correct num- 
ber of WR stars. The largest error is made when using 
single-star model clusters to fit to clusters with multi- 
ple systems. Thus, we choose to fit the observed data 
with simulated clusters containing multiple systems. All 
additional tests described below used simulated clusters 
containing multiple systems and were fit with multiples. 
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3.4.3. IMF Slope 



Fig. 6. — Test-cluster's marginalized ID posterior probability 
density functions (PDF) for a simulated cluster containing only 
single stars that is subsequently fit allowing only single stars left 
or multiple systems right. The input values for the cluster's IMF 
slope, age, mass, and distance are shown as a vertical red line 
in each panel. For a cluster containing only single stars, fitting 
with multiple systems produces extremely biased estimates for the 
cluster age and only slight biases in the IMF slope and cluster 
mass. Without prior information, it is difficult to discern that an 
inaccurate multiplicity model is being used. 

Our bayesian inference methodology shows no system- 
atic biases with respect to the IMF slope as shown in the 
following tests on synthetic clusters. Two clusters were 
simulated with an age of 6 Myr and with IMF slopes of 
a = 2.35 (Salpeter) and a = 0.45 (top-heavy reported by 
IBartko et all (|2010f )). The Salpeter-like cluster has an in- 
put mass of 10,000 M and the top-heavy cluster has an 
input mass of 40,000 M Q in order to produce ~100 non- 
WR stars with Kp<15.5; similar to our observed data 
set. Figure [5] shows the output posterior probability dis- 
tributions for the IMF slope, age, and cluster mass for 
the two simulated clusters. The input IMF slope falls 
well within the 68% confidence interval of the posterior 
probability density function. A handful of similar cluster 
tests were performed with different ages and masses and 
the input and output IMF slopes always agree within 
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Fig. 7. Test-cluster's marginalized ID posterior probability 
density functions (PDF) for a simulated cluster containing multi- 
ple star systems that is subsequently fit allowing only single stars 
left or multiple systems right. The input values for the cluster's 
IMF slope, age, mass, and distance are shown as a red line in 
each panel. For a cluster containing multiple systems, fitting with 
single stars produces extremely biased estimates for nearly all the 
parameters. Without prior information, it is easy to discern that 
an inaccurate multiplicity model is being used since the distance 
tends toward unphysically low values. Multiple systems in clus- 
ters produce bright systems that are brighter than is allowed in 
the model clusters with only single stars at 8 kpc. Thus the only 
way to produce such bright systems is to move the cluster closer 
in distance. 

the 68% confidence region and there appears to be no 
significant bias to either higher or lower IMF slopes. 

4. RESULTS 

The observed Kp magnitudes, their uncertainties, and 
the number of Wolf-Rayet and OB stars were used in 
Bayesian inference to determine the cluster's physical 
properties. Figure |9] shows the resulting ID posterior 
probability distributions for the 4 free parameters in the 
model (cluster age, mass, distance, and IMF slope). The 
model assumes that multiple systems are present with 
the properties described in £|3.1l and Appendix [5] Several 
parameters show correlations including a moderate corre- 
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Fig. 8. — Test-cluster's marginalized ID posterior probability 
density functions (PDF) for a simulated cluster with IMF slope 
= 2.35 (left) and IMF slope = 0.45 (right). The input values for 
the cluster's IMF slope, age, mass, and distance are shown as a 
vertical red lines in each panel. The grey shaded regions represent 
the 68% confidence interval for the output probability distributions 
and all of the input values fall within this interval. The confidence 
intervals are calculated by first finding the peak in the probabil- 
ity distribution and then stepping away from the peak until the 
integrated probability reaches 68%. 

lation between the cluster age and IMF slope fFigureITU]). 
The correlation between the cluster mass and the age or 
IMF slope is a consequence of the age-IMF slope relation- 
ship since, at older ages, the most massive stars have dis- 
sapeared and the cluster mass must increase to match the 
observed numbers of stars brighter than Kp= 15.5. The 
posterior probability distributions are the most accurate 
representation of the results; however, we also present 
"best-fit" values in Table [5] represented by the expecta- 
tion value and 68% and 95% Bayesian confidence inter- 
va 1*3 of the marginalized ID posterior probability density 
functions. The resulting age distribution is multi-modal. 
Therefore, we report solutions from three possible sam- 
ples: (1) ages greater than 3.3 Myr, (2) ages less than or 
equal to 3.3 Myr, and (3) the complete posteriors, The 
age boundaries for the two modal solutions were esti- 
mated based on the minima between the two peaks. We 
report all parameters for these three solutions; however, 
all solutions have similar a, M c ;. (, s , and d distributions. 
We will adopt Solution 1 in the remainder of the paper 
as it contains the most probable solution (with the maxi- 
mum likelihood) and contains the bulk of the probability 
density. Solution 1 is also favored over the smaller age in 
Solution 2 based on more detailed spectroscopic analyses 
that show the Wolf-Rayet star sub-types and the super- 
giant nature of the brightest OB stars favo r ages >4 Myr 
(jPaumard et all 120061: [Martins et al J 120071 ). A compari- 
son of the observed and model KLFs and mass functions 
from the inferred parameters of Solution 1 is shown in 
Figure 111! The complete posteriors for the number of 
predicted Wolf-Rayet and OB stars are also shown in 
Figure [T^] and are compared with the observed number 
of WR and OB stars. The overall fit is good given the 
theoretical uncertainty in the evolution of massive stars 
and the unknown mapping between progenitor mass and 
Wolf-Rayet phases. 

For the IMF slope, recall that a flat prior from 0.5—3.5 
was adopted; but the posterior distribution is strongly 
peaked at 1.7 ± 0.2. The slope is flatter than a typical 
Salpeter slope and slopes of 2.35 or greater are ruled out 

9 Central credible intervals. 
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TABLE 2 

Cluster Properties and Uncertainties 



Cluster 
Properties 

Solution 1: Age>3.3 Myr 

Age (Myr) 

IMF slope (a) 

Integrated Mass (10 3 M Q ) 

Distance (kpc) 
Solution 2: Age<3.3 Myr 

Age (Myr) 

IMF slope (a) 

Integrated Mass (10 3 M Q ) 

Distance (kpc) 
Total 

Age 

IMF slope (a) 

Integrated Mass (10 3 M Q ) 

Distance (kpc) 



Expectation 68% 95% 
Value Interval Interval 



1.2 


3.6 


- 4.8 


3.4 


- 6.1 


1.7 


1.5 


- 1.9 


1.1 


- 2.1 


10.1 


8.1 


11.6 


7.2 


18.5 


7.9 


7.6 


- 8.3 


7.3 


- 8.6 


2.8 


2.6 


- 3.1 


2.1 


- 3.3 


1.8 


1.6 


- 1.9 


1.5 


- 2.1 


8.8 


7.7 


- 9.5 


6.8 


10.5 


7.9 


7.6 


- 8.2 


7.3 


- 8.5 


3.6 


2.8 
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Note. — The integrated masses correspond to that within the field 
of view of the survey (see Paper I for details) between 1—150 Mq . 
The estimated mass for the entire YNC is 2-4 times larger. 
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Fig. 9. — Observed YNC's marginalized ID posterior probability 
density functions (PDF) for the cluster properties. The histograms 
show results from the Multinest bayesian analysis including mul- 
tiple systems. The dashed black lines show the prior probability 
distributions used. The resulting constraints on the IMF slope 
(panel a), the cluster age (panel b) and the cluster mass (panel c) 
are significant compared with the prior PDFs used. The distance 
constraint (panel d) is weak and is largely a reflection of the prior. 

at 3.8<t (99.98%). The slope is inconsistent with the pre- 
vious ly reported top-heavy slope of 0.45 (jBartko et al.l 
120101) at the 3.8a level (99.98%). Our sample of young 
stars includes those within r<0'.'8; however, the results 
are very similar even when these s t ars ar e excluded (Fig- 
ure Hi' as is done in lBartko etHI (|2010[ ). 

The cluster mass is 10,100 M Q (68% interval: 
8,100-11,600 M ) between 1-150 M and is well con- 
strained by the data, given that the posterior distribution 
is significantly different from the prior. The cluster age is 
also significantly different from the prior and shows mul- 
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Fig. 10. — 2D posterior probability density functions (PDF) for 
the observed cluster's properties. The over-plotted contours give 
68%, 95%, and 99% confidence intervals. Weak correlations exist 
between age, mass, and IMF slope. The correlation between the 
cluster mass and the age or IMF slope is a consequence of the age- 
IMF slope relationship since, at older ages, the most massive stars 
have dissapeared and the cluster mass must increase to match the 
observed numbers of stars brighter than Kp= 15.5. 

tiple peaks. The most probable solution falls at 3.9 Myr 
(68% interval: 3.6—4.8 Myr); although the second peak 
at 2.8 Myr is still highly probable. This multi-modal 
behavior is a result of the rapid changes in the number 
of WR stars and the ratio of WR to OB stars in very 
short periods of time (Figure [3j> - Although the distance 
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Fig. 11. — The observed (black) and model (red) Kp-band luminosity function (left) and the resulting mass function (right). The model 
cluster includes multiple systems and has the properties listed in Table l2l (Solution 1). (Left): The observed KLF (black line) and number 
of WR stars (black circle) show good agreement with the model KLF (red line) and model WR stars (red square). The observed Kp 
magnitudes (red) are corrected for differential extinction to a common extinction value of Ak 3 = 2.7. The model KLF is taken as the mean 
and standard deviation (red errorbars) of 100 simulated clusters, all with the same parameters, in order to account for variations due to 
stochastic sampling of the IMF. The imaging completeness curve is applied to each simulated cluster for direct comparison to the observed 
KLF. (Right): The resulting mass function is constructed by converting the observed stars' Kp magnitudes to masses using the best-fit 
isochrone. To capture uncertainties in the stars' brightness, each star's Kp magnitude is randomly sampled 100 times from a Gaussian 
distribution centered on the measured Kp with a width of CTKp- The resulting masses are binned into a mass function and the mean and 
standard deviation of the 100 mass functions are shown as a black line with errorbars. The Wolf-Rayet stars are not included in the mass 
function. 
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Fig. 12. — Comparison of the observed vs. inferred distributions 
for the number of Wolf-Rayet stars (left) and the number of OB 
stars (right). The observed N\yr and Nob are shown as verti- 
cal black lines and the probability distributions from the Bayesian 
analysis arc shown as histograms. 

is also a free parameter, the resulting probability distri- 
bution mainly reflects the prior shown in the dashed line. 
This is reassuring since, in our testing, incorrect model 
assumptions (e.g. fitting single star models to data with 
multiple systems) often led to skewed distance distribu- 
tions, which we do not observe using the multiple systems 
model. 

The cluster mass and the numbers of WR and OB stars 
are based only on the young stellar population within 
the field of view observed for this study, which has in- 
complete azimuthal coverage. This does not represent 
the entire cluster mass of the young nuclear star cluster. 
The total cluster mass can be estimated by extrapolating 
beyond the field of view in this study, in several different 
ways. First, we use the assumption that the cluster is 
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Fig. 13. — Same as Figurc[9] but excluding stars wi th r<0 / /8 to 
better match the sample chosen in jB~a rtko ct al. 20101). The peak 
of the probability distribution for the IMF slope shifts to a =1.6; 
however, it is consistent with the result from the complete sample, 
within the 68% uncertainties. IMF slopes as top-heavy as a = 0.45 
remain inconsistent at 99.89% probability. Ages of 6-7 Myr become 
slightly more probable. 

spherically symmetric and has a surface density profile 
that falls as R~ T with T = 0.9 (paper IT ). This yields a 
surface density weighted coverage of 27% and a coverage- 
corrected total cluster mass of ^37,000 M Q within a pro- 
jected radius of 13'.'65 and above 1 M Q . The above as- 
sumption fails to account for known kinematic structures 
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Fig. 14. — Same as Figure [9] but only allowing single stars in 
the model. The overall fit of the single-star model is poor and 
produces a strong bias away from the measured distance of 8 kpc. 
Furthermore, the resulting evidence for this model is substantially 
lower than for the model that includes multiples. 

in the young population, including a highly inclined disk 
containing as much as 50% of the young stars. Further- 
more, the above mass estimate assumes a constant IMF 
at all radii, which is unlikely given the mass segregation 
that will occur over 4 Myr. As an alternative, we instead 
make an empirically based estimate of the completeness 
by taking the ratio of the number of young stars with Kp 
< 13 within ou r FOV to the same number reported in 
iPaumard et~aTl (|2006f ). which we assume is 100% com- 
plete out to 14". This yields a coverage completeness 
of 71% and a total cluster mass of ^14,000 M Q above 1 
M Q . Thus, we estimate the total cluster mass to be in 
the range of 14,000 - 37,000 M above 1 M . We note 
that the total cluster mass depends on the low-mass cut- 
off as M c i- (mio I 1 Mq ) -1 ' 7 ; and, while we have assumed 
1 M Q , the lowest-mass young stars we observe arc ^8 
M Q . In order to improve the constraints on the total 
cluster mass, more sensitive and azimuthally complete 
spectroscopic surveys at large radii are necessary. 

The above results are derived from modeling clusters 
with multiple systems; however, for completeness, we 
also show the resulting posterior probability distributions 
for a single-star model in Figure Q31 It is clear that the 
fit is poor. The bias to very small distances indicates 
that the brightest OB stars combined with the number 
of WR stars cannot be accurately fit with a single star 
model. Our modeling indicates that multiple star sys- 
tems are necessary to fit the observed stellar population. 
However, we caution that the multiplicity properties were 
fixed to solar neighborhood values that may not be ap- 
plicable to the Young Nuclear Cluster. Both the star 
formation process and subsequent dynamical processing 
may impact multiplicity. 

5. DISCUSSION 

The young nuclear star cluster at the Galactic cen- 
ter, at least along the plane of the well-defined clock- 
wise disk, has an IMF slope in the range a = 1.5 — 1.9 
% confidence interval) for stellar masses above ~10 



M Q . This is significantly flatter than clusters in the 
Milky Way disk and field (a = 2.35). Our measured 
IMF is in good agreement with X-ray observations show- 
ing that X-ray emission from pre-main-sequence stars is 
reduced by a factor of 10 compared with what is ex- 
pected from observations of the Orion Nebular Cluster, 
scaled to the age and dista nce of the Galactic center 
(Navakshi n fc Sunvaevl 120051 ). Extrapolating our mea- 
sured IMF slope down to lower masses produces a factor 
of 4—15 fewer X-ray emitting stars (0.5 — 3.0 M Q ) than 
would be expected for a slope of a = 2.35. 

An extremely top-heavy IMF is inconsistent with our 
observations at high significance, in cont radiction to pre- 
vious observations (|Bartko et al.l 120101 ) . An extensive 
discussion of the possible observational and analysis dif- 
fe rences be tween our work and previous work is presented 
inlPaper J. Here we reiterate a few key differences. The 
IBartko et al.l IpOlOh field of view extends predominantly 
North, out of the plane of the well-defined clockwise disk 
of young stars. They also remove all young stars in the 
central 0'.'8 and outside 12". In the work presented here, 
the spectroscopic window in which young stars are iden- 
tified extends from the center to the ESE along the well 
defined clockwise disk structure. And the entire popu- 
lation of young stars, including those within 0'.'8 of the 
SMBH, is included in the analysis. We assume that all 
the young stars were formed in a single star formation 
event with the same physical conditions giving rise to a 
constant IMF for all the kinematic sub-groups. As dis- 
cussed in Sections [T] and [2j the stars in this innermost 
region were dynamically injected and it still appears to 
be theoretically possible that they originated from the 
most recent starburst that produced the young nuclear 
cluster. Furthermore, restricting our sample of young 
stars to a radius range of 0'/8— 12" does not significantly 
change the KLF shape. Repeating our Bayesian analy- 
sis on this restricted sample produces a similar posterior 
probability distribution for the IMF slopes to that of our 
full sample. 

We find a cluster that is well fit with an instanta- 
neous star formation event at an age of 2.5-5.8 Myr 
with 95% confidence. The most probable age is 3.9 Myr. 
This is smaller than the 6 Myr commonly adopted, al- 
t hough consistent within the uncertainty range reported 
in lPaumard et al.l (|2006r ). The age of the cluster is mainly 
constrained by the number of WR stars and the highest 
luminosity OB stars detected, since this effectively gives 
a main-seque nce turn-off ag e . Ca reful examination of 
Figure 11 in IPaumard et aLl (|2006f) shows that best- fit 
ages, given the different numbers of O and WR stars 
(WN, WC, WNE, WNL), are actually either 4 Myr or 8 
Myr, with 6 Myr being improbable for a Salpctcr IMF. 
The inclusion of fainter OB stars in our spectroscopic 
sample substantially reduces the probability of cluster 
ages as old as 8 Myr. Figure [3] shows the expected num- 
ber of WR and OB stars and their ratio as a function of 
cluster age for different assumed mass functions. While 
our Bayesian analysis captures uncertainties in the age 
due to random sampling of the mass function, it does 
not capture systematic uncertainties in models of stellar 
evolution. The lifetimes and luminosities of WR stars 
are sensitive to metallicity (Fe and a elements), mass- 
loss rates, initial rotation rates, and close-binary mass 
exchange. Increased metallicity, mass-loss rates, rota- 
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tion rates, and binary mass exchange al l tend to increase 
the a ges at which W R stars still exist (lEldridge fc Vinkl 
I2006t iMevnet et al.|[200l ). On the other hand, increased 
metallicity reduces the main-sequence lifetimes. Further- 
more, massive clusters, such as Wcsterlund 1, contain 
WR stars, blue supergiants, yellow hypergiants, and red 
supergiants at an age of 5 Myr; which no stellar evolutio n 
model successfully predicts (e.g. iNegueruela et £l1.II2010T ) . 
Therefore, we caution that ages up to 8 Myr may still 
be reasonable given current uncertainties in massive-star 
evolutionary models. We note that the IMF slope mea- 
surement is dominated by the relative luminosity distri- 
bution of OB main-sequence and supergaint stars, rather 
than the numbers of windy WR stars, and is far less sen- 
sitive to systematic uncertainties in evolutionary models. 
If the young nuclear cluster does have a slightly younger 
age of ^4 Myr, then there are important implications 
for the dynamical history of the cluster. For instance, 
models of in situ star formation with a circular gas disk 
as the initial condition require at least 6 Myr for the 
orbital eccentricities to evolve from circular to the ob- 
served eccentricity d istribution peaked at e = 0.3, assum- 
ing a Salpeter IMF ([Alexander et alJl2007t iCuadra et"afl 
12008t h Revised calculations are necessary to determine 
whether the observed IMF slope and younger age can 
give rise to the observed high eccentricities in such a 
short time. 

Ultimately, a complete spectroscopic and astromctric 
census of the young population in the Young Nuclear 
Cluster is required to accurately measure the recent star 
formation history and determine if the distinct kinematic 
components were formed under different conditions, re- 
sulting in different mass functions. Expanding the sur- 
vey area with current facilities will help clarify whether 
young stars on and off the disk have a different IMF. 
The Bayesian methodology presented here should be ex- 
panded to incorporate additional information from spec- 
troscopy, such as effective temperature and gravity con- 
straints for the OB stars, surface- abundance sub-types 
for the Wolf-Rayet stars, and luminosities/temperatures 
from Wolf -Rayet wind/a t mosph ere modeling such as pre- 
sented in iMartins et al.l ()2007l ) . This additional infor- 
mation may improve age constraints; although uncer- 
tainties in stellar evolution for massive stars will still 
be a major factor. In the future, detection of pre- 
main-sequence turn-on points at Kp=17.5 would more 
definitively constrain the age and star formation history. 
Furthermore, to truly constrain physical models of star 
formation at the Galactic center relative to local star 
formation processes, we will need measurements of the 
IMF shape over a broader mass range, ideally includ- 
ing any turnover suggesting a characteristic mass scale 
(e.g. ~0.5 M Q in Orion). Some theories of star forma- 
tion suggest that in extreme conditions, such as those 
found at the Galactic center, the characteristic mass may 
be significantly higher than in loca l star-forming regions 
(iMorrid H99l iBonnell et all 12004 iKrumholz fc McKed 
12008th Measurements of both the characteristic mass and 
the pre-main-sequence turn-on point should be possible 
with the next generation of large ground-based telescopes 
equipped with adaptive optics-fed integral field spectro- 
graphs. One such instrument is IRIS on the Thirty Me- 
ter Telescope (TMT) , which can accurately spectral- type 
stars in the Galactic center down to Kp=21 in 1 hour 
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Fig. 15. — A model Kp-band luminosity function for the Young 
Nuclear Cluster showing the relative sensitivity of Keck spec- 
troscopy (Kp<15.5) and future TMT+IRIS spectroscopy (Kp<21). 
With TMT+IRIS, the pre-main-sequence turn-on point will be de- 
tectable at Kp~17.5 as well as the full IMF shape down to ~0.4 
Mq . The model cluster shown here has an IMF slope of 1.7 down 
to 0.5 Mq and then tur ns over at 0.5 to a slop e of 1.3 down to 
0.1 Mq as suggested by Wcidncr k Kroupa (2004) for star forma- 
tion in the local neighborhood. The model cluster has a distance 
of 8 kpc and an extinction of Aj{ s = 2.7. 



(|Wright et al.l 12010ft , as shown in Figure [I3J Predicted 
IRIS sensitivities reach 0.4 M stars (Kp=21); however, 
distinguishing the young low-mass objects from the sea 
of much older main-sequence stars may require substan- 
tially higher signal-to-noisc to make precise measure- 
ments of effective temperatures. Such observations will 
be essential as the Young Nuclear Cluster is one of the 
few environments where significant deviations have been 
found from an otherwise near-universal IMF. 

6. CONCLUSIONS 

The Galactic center hosts a young nuclear star cluster 
around the supcrmassivc black hole. We use membership 
probabilities derived from spectroscopy, precise infrared 
photometry, and Bayesian inference methods to deter- 
mine the global properties of the cluster. The best-fit 
age of the cluster is 3.9 Myr, somewhat younger than 
previously reported. The cluster's total mass extrapo- 
lated down to 1 M Q is between 14,000 - 37,000 M after 
correcting for incomplete azimuthal coverage. The best- 
fit IMF slope of 1.7 is flatter than a traditional Salpeter 
IMF; but far steeper than previously claimed in the lit- 
erature. Future spectroscopic observations covering a 
larger field of view and extending to lower masses are 
necessary to understand the impact of mass segregation 
on IMF slope measurements and to determine whether 
the peak in the initial mass function also shows any sig- 
nificant difference from nearby star-forming regions. 
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APPENDIX 
MULTIPLICITY 

The distribution of multiple star systems can be described completely by equations for the multiplicity frequency 
(MF), the companion star frequency (CSF), the mass ratio (q), and the separation distribution. For star clusters at the 
distance of the Galactic center (8 kpc), multiple systems are spatially unresolved and the sepa ration distribution can 
be int egrated over. The multiplicity frequency and companion star frequency are defined as in lReipurth fc Zinneckerl 
(fl99l 

B+T+Q+ 

MF = S + B + T + Q + ... ^ 

CSF = g + 2T + 3Q+ (A2) 
S + B + T + Q + ... v ' 

where S is the number of single stars, B is the number of binaries, T is the number of triples, and Q is the 
number of quadruples. The multiplicity fraction always ranges between and 1; but the companion star frac- 



where CSF>1.5 (jPreibisch et all Il999t IZinnecker fe Yorkd 120071) . Based on observations, the MF an d CSF are 
known to vary with pr i mary star mass, the age of a cluster, and possibly the dens i ty of a cluster ( e.g. iGhez et al 
19931: iPetr et all H998t IMason et all [19981 IDuchene et"aH [l99l IKohler et all l2000t iReipurthl [20001: iPatience et al 
200a IShatskv fc Tokovininl l2002t IKohler et all 120061 : IDuchene et all 120071 : IKohler et all 120081: iMason et all 120091 : 
Raghavan et al.l [20101 iKraus fc Hillenbrand! [20 12). Complete multiplicity surveys of young clusters (<10 Myr) that 
span a large range of primary stella r masses (0.1 — > 10 M^ ), mass ratios, and separations are difficult to condu ct and 
only a few exist in the literature (jLafreniere et al.l 120081 : iKobulnickv fc Frverl 120071 : iKouwenhoven et all I2007D . Star 
formation theories do not yet predict an analytic functional form for the MF, CSF, and q distributions and how they 
vary with pri mary mass: however, simulations a re beginning to produce distributions that are in rough agreement with 
observations (|Batdl2012t iKrumholz et al.ll2012l ). Therefore, we take an empirical approach and compile measurements 
of the MF and CSF as a function of mass from published surveys of young clusters (<10 Myr), which span stellar 
masses of 0.2—17 M and all separations (Table [3]). This compilation is incomplete as we exclude surveys that only 
report multiplicity for masses below 1 M and those without CSF information. The above high-mass limit is set by 
the available multiplicity surveys, and the low-mass cutoff is chosen as a reasonable detection threshold for stars in 
more distant massive young clusters with both current (e.g. Keck) and future (e.g. TMT) spectroscopy. The compiled 
MF and CSF data (Figure IT6|) are fit with a power-law dependence on mass and give the following results: 

MF(m) = Am? A = 0.44 7 = 0.51 (A3) 

CSF{m) = BmP B = 0.50 /3 = 0.45 (A4) 

The best-fit power laws are also shown in Figure [TH For any star cluster older than a few crossing times, the CSF will 
not continue to grow indefinitely at high masses due to the dynamical instability of high-order systems in a clustered 
environment. A hard limit to the mean number of companions is adopted (CSF<3) such that systems with more than 
3 components are allowed, but have a low probability of occurrence. Lastly, some observations sugge st multiplicity 
properties may evolve over time differently for clusters with di fferent masses or stellar densities (e. g. iKroupal 119951 : 
IDuchene et al.llT999t : IKohler et a l. 2006; Mar ks fc Kroupa1l2012[ ). However, the completeness and significance of these 
results is still uncertain ([King et all 120121) and we therefore assume all young clusters (<10 Myr) have the same 
multiplicity properties, including the Young Nuclear Cluster. 
Facilities: Keck: II (NIRC2), Keck: II (OSIRIS) 
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TABLE 3 

Multiplicity Measurements for Young Clusters 



Mass (M Q ) MF CSF 



Reference 



0.175 


0.16 


0.16 


0.390 


0.33 


0.37 


0.915 


0.38 


0.50 


2.14 


0.63 


0.75 


3.4 


0.85 




12.7 


1.00 


1.5 


16.8 


1.00 





Cha I, Lafreniere et al. (2008) 
Cha I, jjafrenic-rcjrt_aL (2008) 
Cha I, l^afrenic-rc ^t_aL (2008) 
Cha I, Jjafrciiicrc^t^aL ^^^) 
Sco OB2. Kouwcnhovcn et al. 
Orion (ONC). Preibisch et ~ 
Cyg OB2. Kobulnickv fe Fryer 




a Mean mass of the sample was estimated from Rizzuto et al. (2011) 
b Mean mass of the sample was estimated from Kiminki et al. (2007) 
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M Q . 



Bate, M. R. 2012, MNRAS, 419, 3115 
Bender, R., et al. 2005, ApJ, 631, 280 

Blum, R. D., Depoy, D. L., & Sellgren, K. 1995, ApJ, 441, 603 
Bonnell, I. A., Vine, S. G., & Bate, M. R. 2004, MNRAS, 349, 735 
Bridges, M., Feroz, F., Hobson, M. P., & Lasenby, A. N. 2009, 

MNRAS, 400, 1075 
Brown, W. R., Geller, M. J., Kenyon, S. J., Kurtz, M. J., & 

Bromley, B. C. 2007, ApJ, 671, 1708 
Carr, J. S., Sellgren, K., & Balachandran, S. C. 2000, ApJ, 530, 

307 

Castelli, F., & Kurucz, R. L. 2003, in IAU Symposium, Vol. 210, 
Modelling of Stellar Atmospheres, cd. N. Piskunov, 
W. W. Weiss, & D. F. Gray, 20P 
Converse, J. M., & Stahler, S. W. 2008, ApJ, 678, 431 
Cuadra, J., Armitage, P. J., & Alexander, R. D. 2008, MNRAS, 
388, L64 

Cunha, K., Sellgren, K., Smith, V. V., et al. 2007, ApJ, 669, 1011 
De Gennaro, S., von Hippel, T., Jefferys, W. H., et al. 2009, ApJ, 
696, 12 

Do, T., Lu, J. R., Ghez, A. M., et al. 2012, submitted to ApJ 
Duchene, G., Bontemps, S., Bouvicr, J., et al. 2007, A&A, 476, 
229 

Duchene, G., Bouvier, J., & Simon, T. 1999, A&A, 343, 831 

Eckart, A., & Genzel, R. 1997, MNRAS, 284, 576 

Eisenhauer, F., et al. 2005, ApJ, 628, 246 

Ekstrom, S., et al. 2012, A&Ain press 

Eldridge, J. J., & Vink, J. S. 2006, A&A, 452, 295 

Espinoza, P., Selman, F. J., & Melnick, J. 2009, A&A, 501, 563 



Feroz, F., & Hobson, M. P. 2008, MNRAS, 384, 449 

Feroz, F., Hobson, M. P., & Bridges, M. 2009, MNRAS, 398, 1601 

Genzel, R., Pichon, C, Eckart, A., Gerhard, O. E., & Ott, T. 

2000, MNRAS, 317, 348 
Genzel, R., et al. 2003, ApJ, 594, 812 
Gerhard, O. 2001, ApJ, 546, L39 

Ghez, A. M., Klein, B. L., Morris, M., & Becklin, E. E. 1998, 

ApJ, 509, 678 
Ghez, A. M., Morris, M., Becklin, E. E., Tanner, A., & 

Kremenek, T. 2000, Nature, 407, 349 
Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 

2005 

Ghez, A. M., Salim, S., Hornstein, S. D., et al. 2005, ApJ, 620, 
744 

Ghez, A. M., et al. 2003, ApJ, 586, L127 
Ghez, A. M., et al. 2008, ApJ, 689, 1044 

Gillessen, S., Eisenhauer, F., Fritz, T. K., et al. 2009a, ApJ, 707, 
L114 

Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009b, ApJ, 692, 
1075 

Goodman, J. 2003, MNRAS, 339, 937 
Goodwin, S. P., & Kroupa, P. 2005, A&A, 439, 565 
Giirkan, M. A., & Rasio, F. A. 2005, ApJ, 628, 236 
Haas, J., Subr, L., & Kroupa, P. 2011, MNRAS, 412, 1905 
Haas, M. R. 2010, PhD thesis, Ph. D. thesis, University of Leiden 
(2010) 



Galactic Center IMF 



17 



Hamann, W.-R., Grafcncr, G., Feldmeier, A., et al. 2008, in 

Astronomical Society of the Pacific Conference Series, Vol. 391, 
Hydrogen-Deficient Stars, cd. A. Werner & T. Ranch, 293 
Hansen, B. M. S., & Milosavljevic, M. 2003, ApJ, 593, L77 
Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377 
Hobbs, A., & Nayakshin, S. 2009, MNRAS, 394, 191 
Hobson, M. P., Jaffe, A. H., Liddle, A. R., Mukeherjee, P., & 
Parkinson, D. 2010, Bayesian Methods in Cosmology, cd. 
Hobson, M. P., Jaffe, A. H., Liddle, A. R., Mukeherjee, P., & 
Parkinson, D. 
Kilbinger, M., et al. 2010, MNRAS, 405, 2381 
Kim, S. S., Figer, D. F., & Morris, M. 2004, ApJ, 607, L123 
Kim, S. S., & Morris, M. 2003, ApJ, 597, 312 
Kiminki, D. C, & Kobulnicky, H. A. 2012, ApJ, 751, 4 
Kiminki, D. C, et al. 2007, ApJ, 664, 1102 

King, R. R., Parker, R. J., Patience, J., & Goodwin, S. P. 2012, 

MNRAS, 421, 2025 
Kobulnicky, H. A., & Fryer, C. L. 2007, ApJ, 670, 747 
Kohler, R., Kunkel, M., Leinert, C, & Zinnecker, H. 2000, A&A, 

356, 541 

Kohler, R., Neuhauser, R., Kramer, S., et al. 2008, A&A, 488, 997 
Kohler, R., Petr-Gotzens, M. G., McCaughrean, M. J., et al. 

2006, A&A, 458, 461 
Kolykhalov, P. I., & Syunyaev, R. A. 1980, Soviet Astronomy 

Letters, 6, 357 

Kouwenhoven, M. B. N., Brown, A. G. A., Goodwin, S. P., 
Portegies Zwart, S. F., & Kaper, L. 2009, A&A, 493, 979 

Kouwenhoven, M. B. N., Brown, A. G. A., Portegies Zwart, S. F., 
& Kaper, L. 2007, A&A, 474, 77 

Krabbe, A., Genzel, R., Drapatz, S., & Rotaciuc, V. 1991, ApJ, 
382, L19 

Krabbe, A., et al. 1995, ApJ, 447, L95 

Kraus, A. L., & Hillenbrand, L. A. 2012, ApJ, 757, 141 

Kroupa, P. 1995, MNRAS, 277, 1491 

Krumholz, M. R., Klein, R. I., & McKee, C. F. 2012, ApJ, 754, 71 
Krumholz, M. R., & McKee, C. F. 2008, Nature, 451, 1082 
Lafrenicre, D., Jayawardhana, R., Brandeker, A., Ahmic, M., & 

van Kerkwijk, M. H. 2008, ApJ, 683, 844 
Laucr, T. R., Faber, S. M., Ajhar, E. A., Grillmair, C. J., & 

Scowen, P. A. 1998, AJ, 116, 2263 
Levin, Y., & Beloborodov, A. M. 2003, ApJ, 590, L33 
Lockmann, U., Baumgardt, H., & Kroupa, P. 2009, MNRAS, 398, 

429 

Lu, J. R., Ghez, A. M., Hornstcin, S. D., Morris, M., & Becklin, 

E. E. 2005, ApJ, 625, L51 
Lu, J. R., Ghez, A. M., Hornstein, S. D., et al. 2009, ApJ, 690, 

1463 

Madigan, A.-M., Hopman, C, & Levin, Y. 2011, ApJ, 738, 99 
Maillard, J. P., Paumard, T., Stolovy, S. R., & Rigaut, F. 2004, 

A&A, 423, 155 
Mafz Apcllaniz, J. 2008, ApJ, 677, 1278 
Mai'z- Apcllaniz, J. 2009, Ap&SS, 324, 95 
Mafz Apellaniz, J., & Ubeda, L. 2005, ApJ, 629, 873 
Marks, M., & Kroupa, P. 2012, A&A, 543, A8 
Martinez, G. D., Minor, Q. E., Bullock, J., et al. 2011, ApJ, 738, 

55 

Martins, F., Genzel, R., Hillier, D. J., et al. 2007, A&A, 468, 233 
Martins, F., Gillessen, S., Eisenhauer, F., et al. 2008, ApJ, 672, 
L119 

Mason, B. D., Hartkopf, W. I., Gies, D. R., Henry, T. J., & 

Helsel, J. W. 2009, AJ, 137, 3358 
Mason, B. D., Henry, T. J., Hartkopf, W. I., ten Brummclaar, T., 

& Soderblom, D. R. 1998, AJ, 116, 2975 
McMillan, S. L. W., & Portegies Zwart, S. F. 2003, ApJ, 596, 314 
Meynet, G., & Maeder, A. 2003, A&A, 404, 975 
Meynet, G, Mowlavi, N., & Maeder, A. 2006, ArXiv Astrophysics 

e-prints, astro-ph/0611261 
Morris, M. 1993, ApJ, 408, 496 
Morris, M., & Serabyn, E. 1996, ARA&A, 34, 645 
Najarro, F., Figer, D. F., Hillier, D. J., & Kudritzki, R. P. 2004, 

ApJ, 611, L105 



Najarro, F., Krabbe, A., Genzel, R., et al. 1997, A&A, 325, 700 

Nayakshin, S. 2006, MNRAS, 372, 143 

Nayakshin, S., & Cuadra, J. 2005, A&A, 437, 437 

Nayakshin, S., Cuadra, J., & Springel, V. 2007, MNRAS, 379, 21 

Nayakshin, S., Dehnen, W., Cuadra, J., & Genzel, R. 2006, 

MNRAS, 366, 1410 
Nayakshin, S., & Sunyaev, R. 2005, MNRAS, 364, L23 
Ncgueruela, I., Clark, J. S., & Ritchie, B. W. 2010, A&A, 516, 

A78 

Nishiyama, S., Tamura, M., Hatano, H., et al. 2009, ApJ, 696, 
1407 

Patience, J., Ghez, A. M., Reid, I. N., & Matthews, K. 2002, AJ, 
123, 1570 

Paumard, T., et al. 2006, ApJ, 643, 1011 

Perets, H. B., Hopman, C, & Alexander, T. 2007, ApJ, 656, 709 
Petr, M. G, Coudc du Foresto, V., Beckwith, S. V. W., Richichi, 

A., & McCaughrean, M. J. 1998, ApJ, 500, 825 
Portegies Zwart, S. F., McMillan, S. L. W., & Gerhard, O. 2003, 

ApJ, 593, 352 

Preibisch, T., Balega, Y., Hofmann, K.-H., Weigelt, G, & 

Zinnecker, H. 1999, ??jnlNew A, 4, 531 
Press, W. H. 1997, in Unsolved Problems in Astrophysics, ed. 

J. N. Bahcall & J. P. Ostriker, 49-60 
Raghavan, D., et al. 2010, ApJS, 190, 1 

Ramirez, S. V., Sellgren, K, Carr, J. S., et al. 2000, ApJ, 537, 205 
Reipurth, B. 2000, AJ, 120, 3177 
Rcipurth, B., & Zinnecker, H. 1993, A&A, 278, 81 
Rizzuto, A. C, Ireland, M. J., & Robertson, J. G. 2011, MNRAS, 
416, 3108 

Sagar, R., & Richtler, T. 1991, A&A, 250, 324 
Salpeter, E. E. 1955, ApJ, 121, 161 
Sanders, R. H. 1998, MNRAS, 294, 35 

Schodel, R., Eckart, A., Iserlohe, C, Genzel, R., & Ott, T. 2005, 

ApJ, 625, Llll 
Schodel, R., Ott, T., Genzel, R., et al. 2003, ApJ, 596, 1015 
Schaller, G., Schaerer, D., Meynet, G, & Maeder, A. 1992, 

A&AS, 96, 269 

Schodel, R., Najarro, F., Muzic, K, & Eckart, A. 2010, A&A, 
511, A18 

Schodel, R., et al. 2002, Nature, 419, 694 

Selman, F., Melnick, J., Bosch, G, & Terlevich, R. 1999, A&A, 
347, 532 

Seth, A. C, Dalcanton, J. J., Hodge, P. W., & Debattista, V. P. 

2006, AJ, 132, 2539 
Shatsky, N., & Tokovinin, A. 2002, A&A, 382, 92 
Shlosman, I., & Begelman, M. C. 1989, ApJ, 341, 685 
Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 
Skilling, J. 2004, in American Institute of Physics Conference 

Series, Vol. 735, American Institute of Physics Conference 

Series, ed. R. Fischer, R. Preuss, & U. V. Toussaint, 395-405 
Tamblyn, P., Rieke, G. H., Hanson, M. M., et al. 1996, ApJ, 456, 

206 

Thies, I., & Kroupa, P. 2007, ApJ, 671, 767 

van Dyk, D. A., Dcgcnnaro, S., Stein, N., Jefferys, W. H., & von 

Hippel, T. 2009, Annals of Applied Statistics, 3, 117 
Veitch, J., & Vecchio, A. 2010, Phys. Rev. D, 81, 062003 
Vollmer, B., & Duschl, W. J. 2001, A&A, 377, 1016 
Wcidncr, C, & Kroupa, P. 2004, MNRAS, 348, 187 
Weidner, C, Kroupa, P., & Maschberger, T. 2009, MNRAS, 393, 
663 

Wright, S. A., Barton, E. J., Larkin, J. E., et al. 2010, in Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 7735, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series 

Yelda, S. 2012, PhD thesis, UCLA 

Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481 



